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Abstract 

The formation of dislocation-free three-dimensional islands during the heteroepitaxial 
growth of lattice-mismatched materials has been observed experimentally for several 
material systems. The equilibrium shape of the islands is governed by the competition 
between the surface energy and the elastic relaxation energy of the islands as compared 
to the uniform strained film. As an exemplification we consider the experimentally in- 
tensively investigated growth of InAs quantum dots on a GaAs(OOl) substrate, deriving 
the equilibrium shape as a function of island volume. For this purpose InAs surface en- 
ergies have been calculated within density-functional theory, and a continuum approach 
has been applied to compute the elastic relaxation energies. 

1. Introduction 

When a materialQis grown on a lattice-mismatched substrate, the uniform strained film 
becomes unstable beyond some critical thickness. Strain relaxation can be achieved 
by the introduction of dislocations. However, there is another important and general 
mechanism of strain relaxation. For Ge/Si(001) Eaglesham and Cerullo[0] observed 
the formation of three-dimensional Ge islands that are dislocation-free. This change 
in surface morphology is driven by the gain in elastic relaxation energy of the islands, 
which overcompensates the energetical cost due to the increase of surface area. Pro- 
vided the interface energy for dislocated islands is sufficiently large this was shown to 
yield the energetically preferred morphology as opposed to both the uniform film and 
dislocated islands for small island size. 101 



*Present address: Physik-Department T30, Technische Universitat Miinchcn, D-85747 Garching, 
Germany. 

^Present address: Department of Physics, Massachusetts Institute of Technology, Cambridge, MA 
02139, USA. 

° Here we consider the situation that the surface energy of the growing material is smaller than 
the surface energy of the substrate, resulting in the formation of a wetting layer and the Stranski- 
Krastanov growth mode. 



Formerly, the roughening of the surface during Stranski-Krastanov growth was con- 
sidered a nuisance for device fabrication. Nowadays, the strain-induced self-assembly 
of small three-dimensional islands is in fact exploited to produce ordered arrays of 
quantum dots.p| A frequently studied examplej^ |, |^, |^, H, |] is InAs/GaAs(001), the 
lattice mismatch amounting to about 7%. Quantum dots have attracted great interest 



due to their zero- dimensional electronic density of states. |Ty, [TT| They may even be of 
technological importance; e.g., buried quantum dots are advantageous for improving 
the device properties of semiconductor lasers. [|1^] 

The dispersion of the sizes of the InAs/GaAs(001) quantum dots is remarkably low.0 
This represents an important and desirable feature of quantum dot growth both in view 
of measurements and device applications, because the delta function shaped density 
of states of every single quantum dot is smeared out by ensemble averaging over the 
size distribution. However, despite its potential importance for the optimization of the 
quality of the quantum dot arrays, the details of the growth mechanism are not yet well 
understood. In particular, the reason behind the narrow size distribution, i.e. even the 
question whether it is due to kinetics or energetics, is still controversial. [|13| , |T^ , p!5[] 

In this article we will focus on the equilibrium shape of coherent three-dimensional InAs 
islands. Knowledge of the equilibrium shape as a function of volume is an essential 
prerequisite for deciding upon the true growth mechanism. Though kinetic effects 
will be important, we expect three-dimensional islands in equilibrium to be observable 
under appropriate experimental conditions: When the concentration of quantum dots 
is low, diffusion of atoms on a single island should be faster than material exchange 
between the islands. Thus shape equilibration should occur on a timescale faster than 
that of Ostwald ripening. Conversely, if already the experimentally observed shape 
deviates from the equilibrium shape, equilibrium thermodynamics will not be adequate 
to describe the island size distribution under the respective growth conditions. 

To calculate the equilibrium shape we have computed InAs surface energies for several 
surface orientations using density-functional theory. The elastic relaxation energy is 
calculated within continuum theory, applying a finite element approach. The total 
energy of the islands, including the strain field in the substrate, is given by the sum 



of these two energy contributions. All the more delicate effects, |jTj] like the strain- 
dependence of the surface energy and island-island interaction are neglected for the 
purpose of predicting island shapes at fixed volume. 



2. Surface Energy 

To determine the InAs surface energies we have carried out total-energy calculations 
within density-functional theory. The local-density approximation is applied to the 
exchange-correlation energy-functional, using Perdew and Zunger's|16|| parameteriza- 



tion of Ceperley and Alder's]!^ data for the correlation energy of the homogeneous 



electron gas. Surfaces are approximated by periodically repeated slabs. The InAs(lll) 
surface, for example, is described by a supercell with a (2x2) surface unit-cell and 
10 atomic layers, the topmost 4 of them being fully relaxed. The atoms in the re- 
maining layers are kept fixed at their bulk positions, using the theoretical bulk lattice 
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constant of 5.98 A. The bottom surfaces of the slabs are saturated with hydrogen- 
hke potentials, [l^ For the Ga-terminated surfaces a Coulomb potential with atomic 
number Z = 1.25, screened by Z electrons, has been used, while a Coulomb potential 
with Z = 0.75 has been taken for the As-terminated surfaces. The saturated surfaces 
are semiconducting without any surface states in the bulk band gap, in this way the 
interaction between both surfaces is made minimal. [|19| For the polar (111) surface of 
GaAs the uncertainty of the total energy due to charge transfer from one side of the 
slab to the opposite side has been estimated to be less than 1.4 meV/A^.|T^ The In 
and As atoms are described by norm-conserving ab initio pseudopotentials[^|, which 
are further transformed into fully separable Kleinman-Bylander ||2 1|| pseudopotentials, 
with the d potential chosen as the local potential. The wavefunctions are expanded 
into plane waves with a kinetic energy < 10 Ry, and the electron density is calculated 
using special k-point sets with the density in reciprocal space being equivalent to 64 
k-points in the whole (100) (1x1) surface Brillouin zone. A generalized version of the 
computer code fhi93cp[^ has been employed. 

To derive a surface energy from a total-energy calculation both the top and the bottom 
surface of the slab have to be equivalent. Such slabs can be constructed for, e.g., the 
(110) and the (100) orientation. The (111) and (111) surfaces of InAs, however, are 
necessarily inequivalent. This is an immediate consequence of the geometry of the zinc- 
blende structure: the bulk lattice can be regarded as a stacking of In-As double layers, 
which are cation and anion terminated towards the top and bottom surface of the slab, 
respectively. Chetty and Martin 1 23 1 solved this problem by introducing an energy 
density. While this energy density itself does not bear any physical significance, they 
showed that integrals over suitably chosen subvolumes of the supercell (e.g., volumes 
bounded by bulk mirror planes) lead to well-defined, physically meaningful quantities. 
In this way the total energy can be divided into a contribution from the upper and the 
lower part of the slab, and after subtracting the respective volume terms we get the 
surface energies of the top and bottom surface separately. This formalism has been 
implemented into the plane-wave code.jlQl 



Similar surface energy calculations for GaAs are described in Ref. [|T^. As the InAs 
surface reconstructions are expected to be equivalent to those found for GaAs, we 
have chosen the same candidates for low-energy surface structures as discussed in that 
reference also for the present calculations. Furthermore, as epitaxial growth most often 
takes place under As-rich conditions, the energies in Tab. |T| refer to As-rich conditions, 
i.e., surfaces which are in equilibrium with bulk As. 

As opposed to the GaAs (110) the As-terminated InAs (110) (1x1) surface does not 
become stable. Independent of the chemical environment, for InAs always the relaxed 
cleavage surface is energetically preferred, displaying the well-known outward rotation 
of the As atom (see e.g. the low-energy electron-diffraction analysis for the (1x1) 



surface performed by Duke et al.\24 



For the (100) orientation our calculation yields the c(4x4) reconstruction as the low- 
est energy surface-structure under As-rich conditions. However, as can be read from 
Tab.0, the energy difference with respect to the a(2x4) reconstruction is so small 
that our calculation is also compatible to the observation of a (2x4) reconstruction. 
Experiment ally |p5| the surface reconstruction has been reported to change from (2x4) 
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orientation 


reconstruction 


surface energy 


meV/A2] 


(110) 


(1x1) relaxed cleavage plane 


41 


(100) 


a(2x4) 


48 


(100) 


c(4x4) 


47 


(111) 


(2x2) In vacancy 


42 


(111) 


(2x2) As trimer 


36 



Table 1: The equilibrium surface reconstructions of InAs under As-rich conditions 
and their surface energies. 



to (4x2) as a function of As chemical potential; in our computations the (4x2) recon- 
struction comes out somewhat too high in energy to give such a transition. Further 
investigations are in progress. 

Both for the (111) and the (111) orientation our predicted equilibrium reconstructions 
are consistent with recent core-level and valence photoemission studies. [p6| , p7[| For 
the (111) the In-vacancy reconstruction is stable independently of the As chemical 
potential. The As-trimer reconstruction, which in case of the GaAs(lll) becomes the 
equilibrium structure under As-rich conditions, is too high in energy for InAs to be 
competitive and thus does not become stable. On the (111) surface again both GaAs 
and InAs display the same As-trimer reconstruction in As-rich environment. 




Figure 1: Equilibrium shape of InAs in an As-rich environment. Surfaces are labelled 
by their Miller indices. 



The equilibrium crystal shape (ECS) of InAs in As-rich environment displayed in Fig.|I| 
has been derived by applying the Wulff construction to the data in Tab. |l[ As we know 
the surface energies only for the {110}, {100}, {111}, and {111} orientations we cannot 
exclude that additional high Miller- index surface orientations may appear on the ECS, 
but the low Miller-index surfaces are expected to be the most prominent ones. As a 
result we find that all four surface orientations co-exist on the ECS, which means that 
they are thermodynamically stable with respect to faceting into each other. This is 
in agreement with the shape of large, and thus presumably fully relaxed, InAs islands 
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grown on a GaAs substrate by metal-organic vapour-phase epitaxy as observed by 
Steimetz et a/.p8| 



3. Elastic Relaxation Energy 

The equilibrium shape of strained three-dimensional islands grown on a lattice- 
mismatched substrate has to be carefully distinguished from the equilibrium crystal 
shape described in the previous section. The optimum shape of a strained island, 
i.e., the shape that corresponds to its lowest total energy, is additionally controlled by 
elastic relaxation. 

We have computed the elastic energy within a continuum theory. The strain field both 
in the island and in the bulk is fully accounted for. For simplicity we have taken the 
experimental second and third order elastic moduli for GaAs[|^ to describe the elastic 
properties of both the island and the substrate. A test calculation with the linear elastic 
constants of InAs and GaAs has shown that this approximation does not affect any 
of our qualitative conclusions drawn with respect to the equilibrium shape. A finite 
element approach is applied to compute the displacement field u(r) and the strain 
tensor e(r): The island and the slab representing the semiinfinite substrate are divided 
into small irregular tetrahedra. The displacement field is defined on the vertices of this 
subdivision, and in between a piecewise linear interpolation of u(r) is used. Thus the 
strain field e(r) calculated from u(r) is constant within each of the small tetrahedra. 
The elastic energy, which is calculated by summing the elastic energy density /(e) 
times the volume over all tetrahedra, is iteratively minimized with respect to u(r). 
Above procedure is repeated for several finenesses of the subdivision, and the energies 
are extrapolated to fineness equal to zero. 

The elastic energy -Emesa of a truncated pyramid with volume V can be approximated 
by a simple analytic expression once the elastic energy of the pyramid with the same 
side faces is known. Let Epy,- denote the elastic energy of this pyramid with volume V. 
/^From the variational property of the elastic energy with respect to the displacement 
field u{r) the elastic energy of the truncated pyramid E' of volume V is estimated by 



with the integral denoting the elastic energy of that part of the pyramid which has 
been sliced off. In the following we are even going to neglect this integral. This still 
yields a good approximation to E' for any not too flat object, because the tops of the 
islands are almost fully relaxed. Finally, we make use of the scaling property of the 
elastic energy, E{V) ~ V, to transform the energy from V back to the island volume 



The quality of this approximation can be judged from Fig. ^. The full lines, given by 
V/V'Epyr, pass through the data points (diamonds and squares, respectively) computed 
for truncated pjnramids with two {111} and two {111} or four {101} faces. 




(1) 



V: 
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4. Equilibrium Island Shape 



To illustrate the basic physical mechanism we are first going to restrict ourselves to a 
very small part of configuration space, i.e., we consider a square based pyramid with 
four {101} facets and compare to the related truncated pyramids which are generated 
by slicing off the top, thus creating a mesa-shaped island with a (001) plane on top, and 
rescaling the lengths such that the volume is kept equal to the volume of the original 
pyramid. 
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Figure 2: Energy gain due to formation of coherent islands (i.e., elastic energy ii^eiast 
plus surface energy Eg^^f minus energy of the uniform strained film Eup) vs. height 
H of the truncated pyramid, normalized to the height i^pyr of the pyramid with the 
same volume. The truncated pyramids are bounded by four {110} facets and one (001) 
facet. The energy curves arc parameterized by the width L of the pyramid, which is 
just a simple measure for the volume. The filled circles denote computed values, while 
the lines represent the analytic approximation discussed in the previous section. 



When a truncated pyramid is created from a pyramid in this way the elastic energy has 
to increase, because material is taken away from the top of the island, where it is already 
almost fully relaxed, and deposited in the remaining mesa where it is still considerably 
strained. On the other hand, given the InAs-ECS from Sect. 2, the surface energy 
decreases. The optimum configuration therefore results from a competition between 
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the elastic relaxation energy and the surface energy. The scaling properties of these two 
quantities, however, are different, the elastic energy being proportional to the volume 
while the surface energy increases with volume like V"^!"^ . Therefore the surface energy 
gains more importance at small volume, while the elastic energy dominates at large 
volume.0 This explains the results displayed in Fig. ^: For small volume the total energy 
as a function of the height of the truncated pyramid has a minimum corresponding to 
some mesa-type shape. When the volume of the island increases, the elastic relaxation 
energy gains more influence and therefore this minimum becomes less pronounced and 
shifts towards the pyramidal geometry. 




1.0 I • ' • ' • ' • 1 

0.40 0.50 0.60 0.70 0.80 

^^JM [meV/A=] 

Figure 3: Elastic energy per volume i^eiast/^ vs. surface energy per volume E^^^ijV 
for InAs islands with volume V = 2.88 x 10^ A^. Squares: square based pyramid with 
four {101} faces and (OOl)-truncated {101}-pyramids. Diamonds: square based pyra- 
mids with two {111} and two {111} faces and (OOl)-truncated pyramids. Triangles 
up: "huts" with two {111} and two {111} faces. Triangles down: square based {101} 
pyramids with {lll}-truncated edges. Dots: islands with four {101}, two {111}, and 
two {ill} faces. Filled symbols denote numerical results, while open circles correspond 
to a simple analytical approximation for (OOl)-truncated "mesa-shaped" islands, as- 
suming that the elastic energy does not change when the (almost fully relaxed) top of 
an island is cut off. Full lines connect islands that are created in this way, varying the 
height of the (001) surface plane. The dashed line is the curve of constant total energy 
-Eeiast + -E'surf that sclccts the equihbrium shape. 

To derive the equilibrium island shape we have to account for the whole variety of 
possible island configurations. To this purpose we have calculated the elastic and 
surface energies for various arbitrarily shaped InAs islands bounded by {101}, {111}, 
and {ill} facets. The results are displayed by filled symbols in Fig.|^. The elastic 
energies of the related (OOl)-truncated "mesa-shaped" islands have been derived by 

"From a practical point of view the volume interval where scaling makes sense is of course limited 
by atomistic effects ioT V ^ and the generation of dislocations at large volume V. 
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means of Eq. (2) as described in Sect. 3, the respective data being denoted by open 
circles in Fig. ^. The optimum shape at given volume corresponds to that point where 
the line of constant total energy touches the manifold of island energies from below. 
Even when the volume is changed, Fig. ^ can nevertheless still be used to derive the 
optimum shape: From the scaling relations we know that the ordinate does not change, 
while the abscissa has to be rescaled according to Esnri/V ~ V~^^^, i.e., only the slope 
of the total energy line decreases when the volume increases. 




Figure 4: The equilibrium shape of a strained coherent InAs island for two different 
volumes, (a) K ~ 8x 10^ A^, (b) V ~ 36x 10^ A^. The orientation of the coordinate 
system is identical to Fig. |I[ 

The equilibrium island shapes (see Fig. ^ are mesa-like hills bounded by {101}, {111}, 
and {ill} facets and a (001) surface on the top, similar to a truncated ECS. Due to 
the different scaling properties of -Eeiast and -Esurf the islands prefer a steeper and more 
pyramidal shape at larger volume while they tend to be flatter at small volume. E.g., 
in Fig. I it can be seen that in comparison to the steeper {111} faces the {101} faces 
are more extended on the small island than on the larger one. Our equilibrium island 
shapes differ from the {101}-pyramids grown experimentally by Ruvimov et a/.[0 for 
yet unknown reason. They are, however, similar to the shapes of InP islands on GalnP 
observed by Georgsson et a/.[^ 

The equilibrium island shape evolves continuously with respect to the volume. There- 
fore the total energy of the equilibrium islands does not simply vary with volume 
like a + bV~^^^ as would follow from simple scaling. Instead, when we represent the 
low-energy envelope of the data in Fig. || by y = a + bx~"', with x = Esuvf/V^^^ and 
y = Ee\a.st/V, and a, b, and 7 independent of volume, b and 7 positive, we get the total 
energy per volume from the Legendre-transform of y{x): 

E,jV = a + b'V-''/'^^+'\ (3) 

with some positive constant b', i.e., the volume-exponent becomes smaller. Of course 
our approach which, e.g., disregards surface stress effects, still leads to Ostwald ripening 
by means of its construction. To discuss the possibility of an optimum island size one 
has to consider further elastic interactions. 
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